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Abstract 

We discuss a family of planar (two-dimensional) systems 

with the following phase strucure: a Fermi liquid, which goes by 

a second ordertransition (with non classical exponent even in 

mean-field) to an intermediate, inhomogeneous state (with 

nonstandard ordering momentum) , which in turn goes by a first 

order transition to a state with canonical order parameter. We 

analyze two examples: (i) a superconductor in a parallel 

magnetic field (which was discussed independently by 

Bulaevskii)for which the inhomogeneous state is obtained 

for 1.86T C < B < 1.86\/2T C where T c is the critical temperature (in Kelvin) 

of the superconductor without a field and B is measured in 
Tesla, and (ii) spinless (or, as is explained, spin polarized) 
fermions near half-filling where a similar, sizeable window (which 
grows in size with anisotropy) exists for the intermediate CDW 
phase at an ordering momentum different from (tt, vr). We 
discuss the experimental conditions for realizing and observing 
these phases and the Renormalization Group approach to the 
transitions. 
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I. INTRODUCTION 



Zero temperature quantum critical phenomena have been a subject of intense theoretical 
and experimental investigation in the recent years. Of particular interest to us has been 
the Renormalization Group (RG) description of the transitions. In our pursuit of the RG 
description of two-dimensional fermionic systems, we came across a family of systems which 
exhibit similar phase structure and similar phase transitions separating the phases in mean- 
field theory. Since this phase structure seems generic and experimentally accessible today, 
we have chosen to make it the focus of this paper, only briefly discussing the RG program 
that led us to them in the first place. 

The systems in question exhibit three phases as a control parameter / is varied. Since 
the physics depends on the ratio W/I, where W is the interaction strength, we can either 
vary / at fixed W (which is what happens in experiment), or vary W at fixed I in a theoret- 
ical discussion. The phases of the system are a weak coupling (large /) Fermi liquid phase, 
and a strong coupling phase with conventional order, separated by an intermediate, inho- 
mogeneous phase with nonstandard ordering momentum. The intermediate phase, which 
is the most interesting, is accessible over a broad window in I (in clean systems) thanks 
to two-dimensionality. Indeed, in anisotropic systems that are more one-dimensional, the 
window is broader. 

We will illustrate our theme using two concrete examples: a system of electrons with 
an attractive coupling between up and down spins and a system of spinless fermions with 
nearest neighbor repulsion. (We hasten to add that we will indicate how these "spinless 
fermions" are to be realized experimentally by polarizing ordinary electrons in a manner to 
be specified.) 

Consider first the system of free electrons - by which we always mean Landau's quasipar- 
ticles - with a circular Fermi surface (FS). They are unstable to the BCS state at arbitrarily 
small attraction W. The instability may be traced back to the time-reversal symmetry of 
the spectrum and hence of the FS: E(K) = E(—K). 
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Consider next a system of free spinless fermions, at half-filling, with nearest neighbor 
hopping on a rectangular lattice. This system is unstable to arbitrarily weak nearest neighbor 
repulsion W which drives it to a Charge Density Wave (CDW) state at momentum Q = 
(7r,7r). The immediate instability is again due to a special symmetry of the FS: due to 
particle-hole symmetry and nesting, if K lies on the Fermi surface, so does K + Q. 

Both 

instabilities have a nice RG description. It was shown in Reference Q that in an RG 
scheme in which modes around the Fermi surface are symmetrically eliminated in thin slices, 
the Fermi liquid appears as a fixed point and the above mentioned interactions appears as 
a marginally relevant terms. The existence of a steady flow of the coupling in the simple- 
minded mode elimination scheme was also seen to result from the special symmetry of the 
Fermi surface. 

We started by asking if such an RG description could be extended to transitions at 
non-zero coupling . A necessary prelude to such an enterprise was finding systems which 
exhibited such transitions. It was at that stage that we ran into the systems discussed here. 

Consider first the superconductor. To move the transition to a nonzero, but small cou- 
pling, (which will allow a perturbative RG analysis ), we must destroy the time- reversal 
symmetry of its FS by a small amount. To this end, let us apply a magnetic field B parallel 
to the plane. This parallel field has no effect on the orbital motion but causes a Zeeman 
splitting of size 2fiB = 21, where /i is the electron's magnetic moment. We expect that 
at very small fields, the superconductor will be stable, while for very large fields the Fermi 
liquid will be stable. 

Our strategy for showing that there exists an intermediate inhomogeneous phase is as 
follows. We start at the BCS end and ask when the BCS state yields to the Fermi liquid 
as the field is raised. Then we go the Fermi liquid end at large / and ask what happens to 
it as the field is lowered. We find that it first becomes unstable to inhomogeneous pairing 
before the BCS state beats it in energy. There is clearly a window when neither the BCS 
nor the Fermi liquid can be the ground state. Here the inhomogeneous state prevails and it 
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is to this state, and not the Fermi liquid that the BCS state succumbs as the field is raised. 

Here are some details on the implementation of our strategy. Consider the response of 
the BCS paired state to a small magnetic field B. Since the spins form a singlet they ignore 
the field altogether. In the meantime, as B increases, the Fermi liquid energy decreases 
because it is getting polarized. Beyond the Chandrasekhar-Clogston [F2J limit Ip = A /y/2 
(where A is the BCS order parameter in the absence of the field), the Fermi liquid has lower 
energy than the BCS state and one may expect a first order transition (hence subscript F 
on Ip ) to the polarized Fermi liquid (PFL). But what happens is that by this point a state 
with inhomogeneous order has beaten the Fermi liquid, and the first order transition is really 
to this state. 

To see all this, so let us go to the weak coupling end (large /) where it is natural to first 
find the FS in the presence of the field and then turn on the attraction. The radii of the up 
and down FS's now differ by q = 21 /vf, where vf is the Fermi velocity, hereafter set equal 
to unity. The PFL is stable to pairing fluctuations, at least for small coupling or large /. 
If we perform the sum over Cooper bubbles and evaluate the pairing susceptibility as I is 
lowered, we find two interesting features: 

• The PFL susceptibility first diverges at I g = Aq, which is \/2 times bigger than the If 
at which the BCS energy rises above that of the PFL. 

• This pairing is at momentum q = 21 . 

Thus we have a window 

A /V2 = I F < I < A = I q 

when the BCS state has a higher energy that the PFL, which itself is unstable to the in- 
homogeneous state. To repeat, as I is increased, the BCS state will yield, by a first order 
transition, not to the PFL, but to an inhomogeneous state with wave number q ~ 21. The 
phase diagram is shown in Figure la. Actually the inhomogeneous state will prevail even 
slightly below the lower limit of the above window, since at the first order transition the 
BCS state has to beat, not the PFL, but the inhomogeneous state 
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which has a lower energy. (For completeness we mention that at a still lower field, 
Jo = A /2 the q = susceptibility of the Fermi liquid diverges. This is however not very 
significant since the corresponding instability is preempted by the other transitions. ) 

Now for the spinless fermions. Here we destroy the perfect nesting of the FS by introduc- 
ing a chemical potential I. (This possibility was raised in Reference |T[ along with another 
option: adding a second neighbor hopping term, which we will consider only briefly.) We 
establish the window as before. Let us start at strong coupling (small I), and consider a 
CDW state at half-filling, with an order parameter Ao- Now imagine modifying the chemical 
potential. The system does not respond to a change in the chemical potential which is mov- 
ing in the gap, exactly the way the superconductor did not respond to the applied field. Once 
again a mean-field calculation shows that the (jr, it) CDW and the FL have equal energies 
at Ip = A /-\/2, where A is the condensate without the chemical potential term. . At the 
weak coupling end, the FS has shrunk, as shown in Figure 2, and becomes stable against 
particle-hole pairs condensing at Q — (it, it) or any other momentum. An RPA calculation 
shows that as I is lowered, the FL becomes unstable to a CDW at a momentum different 
from (tt, 7r) at I p rs rA , where r > 1 is an anisotropy parameter which measures the ratio 
of the hopping in the y and x directions. Thus we have a window Ao/y/2 < I ~ rA when 
an inhomogeneous state must prevail. This state is inhomogeneous in the sense that it is 
characterized by a momentum different from the canonical value, which in this case happens 
to be (vr,7r). The phase diagram is shown in Figure lb. (The significance of the points I q , 
and Iq which do not correspond to transitions, will be explained in section III.) 

In both the BCS and CDW cases the mean-field second order transition from the FL to 
the inhomogeneous state was found to have non-classical exponents. For example (3 = 2 for 
the superconductor. This nonclassical answer in a mean-field calculation may be ascribed 
to the fact that the mode that becomes unstable has a momentum which is a singular point 
of the susceptibility. This singularity and the existence of a fairly large window for the 
inhomogeneous phase are both a result of the two- dimensionality of the system. 

The plan of the paper is as follows: In the next section we will discuss the superconductor 
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and relate it to earlier work. In trying to track down existing literature in the case of the 
superconductor, we were invariably directed to the works of Larkin and Ovchinnikov, and 
Fulde and Ferrel ||, and we shall briefly discuss these. However, upon digging further, we 
found some work of Bulaevskii [Q to be far more relevant, as will be discussed. In section 
III we will perform a similar analysis of the spinless fermion system and describe how it may 
be experimentally realized. Section IV is devoted to the status report of the RG program. 
We end with conclusions and avenues for future research in section V. 



II. THE PLANAR SUPERCONDUCTOR 

We will begin with noninteracting electrons lying within an annulus of thickness 2A 
centered around the Fermi circle of radius Kp. We will construct a mean- field hamiltonian 
that pairs the up and down electrons whose momenta lie within this band. 

The energy of an up electron is 

E + (K) = I + ( R2 ~J {2p ) = I + v F k = I + k (1) 
where the Fermi velocity Vp has been set equal to unity, and 

k = K — Kp. (2) 
Terms higher order in k are dropped. 

— * 

The energy of its partner, a down spin electron of momentum — K + q, is 

E_(-K + q) = -I + k-qcos6 (3) 

where 9 is the angle between K and q. 

(It is understood that — K + q is also constrained to lie in the same annulus. This means 
that certain electrons will not have partners for pairing in a state of momentum q. They 
must be handled correctly.) 

A useful combination is 
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E + (K)-E_(-K + q) 



I+(q/2)cos9 = z(9). 



(4) 



Let us begin with the mean-field hamiltonian 
AA* 



H 



W 



//V p *2tv 

dk d6 \c ] + (K)c + (K)(I + k) + cl(K)c_(K)(-I + k) 
-A J L 

-6(A - ||i? - q\ - K F \) \Acl(-K + q)cl(K) + A*c_(i?)c + (-i? + g*) 



(5) 



where the condensate was taken to have the form A(x) = A exp [iq ■ r\ , the step function 
0(A — 1 1 if — q\ — Kp\) ensures that the other partner in the pair also lies in the annulus, 
and the subscripts refer to the spin. The reader should not worry about factors of 2n since 
all key results will be given in term of observables. 

By going to Bogolubov operators d± and performing the usual transformation that kill 
the "bad" terms, we obtain 

W 

+ [ dk [ n d6 0(A- \\K-ql -K F \) 

J-A JO 



k- I -qcos9 + dld + I + (q/2) cos 6 + A 2 + (k - (q/2) cos 6) 
+ {dld_ - 1) [-/ - (q/2) cos6 + \j A? + [k - (q/2) cos#) 2 



+ ^--4A 9 . 



(6) 



The last two terms come from the unpaired electrons. The ground state is found by filling 
all negative energy states. (If we set A = above, the ground state energy should reduce 
to — 27r(A 2 + J 2 ), the energy of the PFL. ) If we subtract off the energy of the PFL, we find 
the relative energy 



A 2 n 2 2A 
E(A) - E(0) = 2vr/ 2 + — + 7r(±- - A 2 ) - 2ttA 2 In 

W 4 A 



r2vr 

+ / dO 



Jo 



e(z-A) 



a2, z + \/ z 2 - A 2 

-zVz 2 - A 2 + A In 



+ z 



(7) 



z(6) = I + j-cos6 



where A and A 2 stand for |A| and |A| 2 . 



Most of our results are based on an analysis of this equation. 

The first term comes from subtracting off the polarization energy of the PFL. As for the 
9 integral let us note that if A > I + q/2, it does not contribute. It is clear that nonzero q is 
a liability in this region and the minimum energy is readily found to occur at A given by 

exactly as if B were never turned on. This is just the manifestation of the fact that the BCS 
system ignores the field. The condensation energy (relative to the PFL energy) is found to 
be 

E(A) - E(Q) = vr(2/ 2 - Aq) (9) 

which leads to the result that I F = A /\/2 is the field at which the BCS system yields to 
the Fermi liquid. (Since A = y/2Ip > Ip the answer is consistent with the assumptions 
made in deriving it.) 

Now for the weak coupling (or large J) Fermi liquid end. Either by summing all the 
Cooper bubbles, or by taking the A 2 derivative of the energy at A = 0, we find the formulas 
for (In taking this derivative it is important to note that if I > (q/2), z{9) is always 
positive and only the first term in the 9 integral contributes, whereas if {q/2) > I, the second 
contributes for A = 0.) The result is 

x^^w-jy^m (10) 

1 , 2A 
2tt In . 

(1 I + y/P- (q/2) 2 

2A 



1 2A . . . 

— -27rln (?/2) </ (11) 



1 2A 

7I7- 27rln 777^ (9/ 2 ) > L ( 12 ) 



(13) 

Note that at each value of I (the applied field) , the softest mode is the one with (q/2) = I 
and that this is a singular point of x- This mode becomes unstable when / = I q , where 



If we compare this equation to the gap equation in the absence of field 



1 




(15) 



W 



we see that 



A, 



o 



(16) 



Thus we have a window 



=L < I < A: 







(17) 



wherein the neither the Fermi liquid, nor the homogeneous superconductor is the ground 
state. (Actually the lower limit should be lower since the BCS state has to beat not the 
PFL but the inhomogeneous state. However the condensation energy of the inhomogeneous 
state is quite small and so is the change in the lower limit.) 

If we put in the numbers we find that if B is in Tesla and T c in Kelvin, the new phase 
should be stable in the interval 1.86T C < B < 1.86y/2T c . We have used here the free electron 
g factor. In practice g could be much higher, and the values of B correspondingly lower. 

Figure 3 shows a plot of the energy in Eqn.(^) versus A. The parameters A = 100, / = 
.834, W = .031 are chosen so as to bring the BCS energy (measured relative to the PFL) 
to zero, i.e., we have chosen / = A /a/2. Plotted on the same graph is the energy of the 
inhomogeneous state at (q/2) = .923. Note that this (q/2) is not equal to /. In other words, 
although (q/2) = I becomes unstable first, once the order sets in, a larger (q/2) does better 
and the graph shows the energy for the best (q/2). If we lower I a bit, the BCS state will 
dip down very quickly to below the inhomogeneous state. Thus we may take I = Aq/^/2 to 
be essentially the field for the first order transition to the inhomogeneous state. If I is raised 
further, the BCS graph will move upwards and so will the inhomogeneous state's energy 
at the minimum. As the minimum moves towards the origin, the optimal (q/2) will move 
towards the current value of /. Finally for / > Ao the PFL will become stable to all pairing 
fluctuations. The transition to the PFL is clearly second order in mean-field theory. 
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From the energy formula it is easy to compute the specific heat by finding the density 
of states, and the magnetization by taking the I derivative. (Since (q/2) and A are slaved 
to I one may ask how they are to be handled during the derivation. In general there will 
be implicit and explicit changes in E to to a change of /, however at a point where dE/dq 
and dE/dA vanish (as they do for us) only the explicit derivative matters.) The result is, 
in terms of the corresponding quantities for the PFL, 

wt-r^iC deV ^ Wz - A) - e( - z - A)] (18) 

C 1 j 2 * m [*H z - A) - z6(-z - A)] 



Cpfl 2n JO \J z 1 — A 2 

where A and (q/2) are at their optimal values. Putting in the numbers we find that these 
ratios drop from the value of unity at the PFL end, down to .92 and .79 respectively by 
the time we reach the situation shown in Figure 3 (which is more or less the end of the 
inhomogeneous phase). 

Let us next compute the exponent j3 which characterizes the onset of A. If we could 
write the energy near the transition as 

E = -tA 2 +uA 4 + ... (20) 

where t = A — /, we would find as usual that A ~ t 1 ! 2 . If however we analyze our energy 
function we find that 

E = -tA 2 + uA b ' 2 + ... (21) 

leading to (5 — 2. In other words, the energy function does not have an analytic expansion 
in A near the transition, which in turn can be traced to the fact that the optimal (q/2) at 
the transition is a singular point of x^ 1 ■ Thus even though we are doing mean-field theory, 
a nonclassical exponent emerges. 

We could go on and compute the electromagnetic response functions but do not, for the 
following reasons. Our analysis has shown that in the window 
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an inhomogeneous state characterized by momentum of magnitude q (which is of the order 
/) has lower energy than either the PFL or the BCS state. However our analysis does not 
choose a particular direction for q. Thus, as in problems involving, say crystallization, 
all modes in this circle are degenerate. In this case the low energy physics of the condensate 
is very much like that of a Fermi system, i.e., centered around a circle and not the origin. 
Recently Hohenberg and Swift have shown how the shell elimination devised in Ref. 
for fermions can be applied to this problem. For the present simply note that, instead of 
the plane wave we considered so far, a standing wave of the same wavelength might have a 
lower energy. The standard trick for deciding what exactly happens is to expand the energy 
functional in a power series in A(g), going out to fourth order. The fourth order vertex 
(which is the particle-particle bubble with two more insertions it it) will determine which 
combination of modes is best the near transition. Unfortunately we cannot do this here 
since an analytical expansion in A does not exist. While we may expect that (3 = 2 may be 
insensitive to the actual form of the order parameter (plane wave versus cosine) , features like 
specific heat will be very sensitive to the true ground state. For example if A ~ cosgr, bands 
will be formed and the density of states can go up above the PFL value. The only thing that 
seems certain is that the specific heat will not be activated. Likewise the electromagnetic 
response and Meissner effect (for a small perpendicular field) will depend on the knowledge 
of the true ground state. This is a problem we have not solved. All we can say is that if the 
answer is a cosine, there will be an oscillation in spin density at the same wavelength. 

We now consider the relation of our work to that of Larkin -Ovchinnikov and Fulde-Ferrel 
|| These authors considered a three dimensional superconductor with ferromagnetic 
impurities. The spins of these impurities coupled to the electrons via the exchange interaction 
so that they could be represented on the average by some external field that coupled only 
to spin. These authors too predict an inhomogeneous phase. We now list the differences 
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between our work and theirs. 

• These authors consider a system coupled to magnetic impurities while we need a clean 
system. 

• They evade the orbital complication of the magnetic field (even though they work 
in three dimensions) because the field really represents the exchange interaction with 
impurities. We evade it by considering a planar system with an external parallel field. 

• Most importantly, they expect the inhomogeneous phase over a very narrow window 
.707Ao < / < .755Ao in three dimensions whereas we expect it over a broad window 
.707Ao < / < Aq in two dimensions, with I as the controllable external field. The 
window will be even wider in anisotropic systems (e.g. elliptical FS) because they are 
more one dimensional and will have roughly parallel segments of the FS over longer 
intervals. 

• Due to the fact that the optimal q is a singular value of \~ l we find nonclassical 
exponents {(3 = 2) in the mean-field level, while they find classical exponents due to 
the fact that the optimal q, for small A is q = 1.21, which is a nonsingular point of 
the three dimensional susceptibility. This allowed Larkin and Ovchinnikov to study 
mode coupling and to decide that the cosine order parameter is better than the plane 
wave. Till we can do a similar thing here it may be a reasonable to assume that the 
same thing happens in d = 2 also. In this case we can expect as they do, that the spin 
density will will oscillate with period q and that the specific heat will be greater than 
for the PFL. 

Consider next the work of Bulaevskii ||. He looks at layers in a magnetic field and 
works out critical fields for any tilt and at finite temperature. His formula for zero tilt at 
zero temperature gives the window we quoted. The exponent (3 = 2 is reflected in the T 
dependence of his results. Like us, he does not consider spin orbit coupling and for this 
reason could not explain some of the experiments || done at that time, as pointed out by 
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Klemm, Luther and Beasley || , who did include this effect. It appears that the time is ripe 
to see this inhomogeneous phase. The conditions are that the material be as two-dimensional 
as possible, clean, have low Z (to minimize spin-orbit effects) and have a low T c so that the 
required fields are are not too large. It will be interesting, not only in its own right, but 
also as a part of the family of systems discussed here. Of course at finite T, there will be no 



ordered phase and a Kosterlitz-Thouless phase with algebraic order will take its place [|10 



The phase structure will still be visible at small T and the T = critical point will control 
the finite T physics in the quantum critical region as in other problems |l 1 1 (12 . 



III. THE SPINLESS FERMION SYSTEM WITH CDW ORDER 

We begin by reassuring our readers that what we mean by a spinless fermion can be 
experimentally realized in at least one way. In the idealized calculation we consider a 
system of spinless fermions that fills up a nested Fermi sea and ask what happens as we 
tamper with the nesting by changing the chemical potential. Assuming that all the action 
is centered around the immediate vicinity of the Fermi surface, i.e., that the interaction and 
the change in chemical potential /, are small compared to the bandwidth, we can duplicate 
the same physics near the Fermi surface as follows. We take regular electrons and fill them 
up to, say, 45% (instead of half) of the Brillouin zone. Now we apply a magnetic field in 
the plane. This does nothing to the orbital motion but splits the spin up and down Fermi 
surfaces. For some choice of field, the bigger of the surfaces will begin to nest, while the 
other would have moved away in the opposite direction and hopefully will not do anything 
interesting. The polarized electrons with a nested or nearly nested Fermi surface will be 
the spinless fermions of the theoretical calculation. Of course a CDW in this case will also 
imply a spin density wave or SDW. 

As in Reference |jj we will consider free fermions with the dispersion relation which 
comes from nearest neighbor hopping: 

E(x, y) = — cosx — r cosy (22) 
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where x and y stand for K x and K y . We choose anisotropic hopping since the generic 
physics we are talking about can be obscured by the van Hove singularities that lie on the 
Fermi surface when r — 1. We will work with r > 1. This energy relation satisfies the 
condition 

E(x + TT,y + Tr) = -E(x,y) (23) 

and the vector Q = (jr, it) connects points on the two branches of the Fermi surface, 
labeled by a = ±1 as shown in Figure 2. At half-filling the Fermi surface is defined by 
E = or 



y — ± cos 1 



cos a; 



a cos 1 



cos a; 



(24) 



We will shortly use a variable e that measures the energy from the FS of the free particles. 
In our problem where the FS is at zero energy, e = E. 
Imagine turning on a chemical potential so that 

E(x, y) = — cos x — r cos y + 1 (25) 

The Fermi surface of the free fermions will now shrink to a smaller size and no longer nest. 
Of relevance is the combination 

E(x + 7T + q,y + 7r + p) -E(x,y) _ , , . . 
2 = z{x ) y ) q,p) (26) 

= / — (q/2) smx — r(p/2) siny (27) 



= I — (q/2) smx — (p/2) aV r 2 — cos 2 x (28) 
= z a . (29) 

In the above formula, a = ±1 is the branch index, and, as everywhere else in the paper, 
terms of quadratic order or higher in the small quantities I,p, q will be ignored. 

Once again we limit ourselves to a band of energy A on either side of the free electron 
FS. To this end we change variables from (x, y) to (x, e), where e measures the energy from 
the FS. (See Reference for details. ) As a result 
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J dx dy —> J dxde.J(r, x) 



(30) 



where the Jacobian is 



J(r,x) = 



\Jr 2 — 



cos 2 x 



(31) 



on the free electron FS. We will ignore its deviation off the FS. It is useful to define 

_ 1 r 27r 



— 1 r Z7r 
J = — / J(x)dx. 
2tt Jo y ' 



(32) 



The mean-field hamiltonian is 

A 2 f A f 2lT r t t 

H — — + / de \ dx.J(r,x) c+(e, x)c+(e, x) + cL(e, x))c-(e, x) 

W J -A JO 

/A p2tt 
de / dxJir, x)Q(\e(x + ir + q, y + it + p)\ < A) 
-A JO 

A(c\.{K)c-(K + Q + iq + jp) + A*cl(K + Q + iq + ]p)c + {K) 



(33) 



where now the subscripts ± on the operators refer to the branches of the FS. 

To exhibit the intermediate phase, let us assemble the two ingredients: (i) The energy 
difference between the Fermi liquid and the (it, it) CDW state and (ii) The inverse Fermi 
liquid susceptibility (in RPA) to CDW formation at momentum (x + it + q, y + tt + p). The 
space between the zeros of these gives a window in I for the inhomogeneous state. 

A standard mean-field calculation, as in the superconducting problem, shows that for 
A > I, the energy difference between the (tt, tt) CDW state and Fermi liquid is 



A 2 /-27T A 2 
E(A) — E(0) = — + I dxJ(r, x) [I 2 - — - A 2 



In 



2A, 



(34) 



This energy difference has a minimum at 

A = 2A exp 



1 



(35) 



2ttJW. 

Feeding this back into the energy formula we find that the / at which the (tt, tt) CDW state 
loses out to the Fermi liquid is 

r A 



V2' 
16 



(36) 



(Note that once again the answer is consistent with the assumption A > I that went 
into its derivation.) 

By summing the usual particle hole bubbles, the inverse susceptibility of the Fermi liquid 
is found to be 

A 



1 1 r 2n A 
~\<1,P) = ttt - / dxJ(r,x)\n — 

W 2 „ JO \Za\ 



X 



where z a is the value of z on branch a, as defined in Eqn.(|29|). This gives 

1 r 2n A 
X _1 (g,p) = — -/ dxJ(r, x) In 

W Jo ^(J - (q/2) sinx) 2 - (p/2) 2 (r 2 - cos 2 x) 



(37) 
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First note that the uniform CDW (i.e., (tt, it)) mode becomes singular when I = Iq where 

(39) 



i- = 27rJln^ 
W I 



i.e., at 



Aexp 



1 



Ao 
2 ' 



lixJW. 

Next we consider just nonzero q. The instability occurs when 

A 



(40) 



- = J(r, x)dx In — 



(41) 



(q/2) sinx| 

Let us consider this equation at large r where the analysis is easier. Now J ~ 1/r and 



1 2i, 2A 
— = — In -j= 

W r I + yJP 

2tt 2A 
= — In 



(?/2) 2 
q/2 > I. 



q/2<I 



(42) 
(43) 



r (q/2) 

The optimal value is q/2 = I, and we find upon inverting the above that the I value at 
the onset of the instability is 

1 



2Aexp 



2nW 



2 Aexp 



2-nJW 



An 



(44) 



At r not too large, but also not too close to unity (so that the singularity of the Jacobian 
does not dominate the integral), similar results hold. 
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So far we have found the same window as in the case of the superconductor, namely 

A 



V2 



< I < Ao 



(45) 



But the window gets wider if we look at nonzero p. (This is also clear from Figure 2.) 
Now we have 



1 r Z7r A 
— = J(r,x)dx\n = 
W Jo jp_ ( p /2) 2 (r 2 -cos 2 :r) 

= / J(r, x)dx In ■ 
J o 



(46) 



Jp _ (rp/2) 2 + (p 2 /8) + (p 2 /8) cos2a;) 

For large r, this integral is seen to be a maximum when p = 21 /r where its value is 
(2vr/r) ln^ so that 



I p = rA . 



(47) 



The following table shows the deviation from this result for smaller r. In all these cases 
the optimal p was found to be 27/r. Note that the large r results are not bad for r as low 
as 1.1. 



r 


1.1 


1.5 


2.5 


3.5 


rAo 


.83 


.93 


.98 


.99 



The final phase diagram is as shown in Figure lb. The main point is that as the anisotropy 
grows, so does the window, since the system is becoming increasingly one dimensional and 
the opposite faces of the Fermi surface are becoming globally parallel. (In the limit of parallel 
faces there exists a perfect nesting vector for any I.) 

We will not consider the general case of a nonzero q and p since we have already seen a 
substantial window for observing the inhomogeneous phase. We also do not try to compute 
its energy as a function of A (as we succeeded in doing for the superconductor) because of 
the following problem. Recall that in the superconductor problem we coupled K to — K + q. 
That in turn must be coupled to — (— K + q) + q = K, which means only two modes are 
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to be coupled. In the present case we coupled K to K + (ir, ir) + q, which in turn must be 
coupled to K + 2q and so on. We did perform a calculation keeping just the first coupling 
to exhibit a variational wave function which could beat the Fermi liquid and the (ir, ir) state 
inside the window. 

Due to the fact that the optimal momentum p at the transition is a singular point of x, 
we can once again expect nonclassical exponents. We guess that (5 will be larger than the 
classical value, probably equal to 2 once again. 

A. Experimental realization 

Let us now turn to the experimental realization of this system. We start with a system 
of real electrons not too close to half-filling, say at 45% filling. At zero temperature, this 
will probably be a BCS superconductor. Next we apply a parallel magnetic field. This is 
the problem we just studied. The BCS state will eventually pass via the inhomogeneous 
superconducting phase to the polarized Fermi liquid as the field rips apart the Fermi surfaces 
of the up and own spins. (Here we must choose a material with a large g factor for the 
electron.) Focus on the two Fermi surfaces. One shrinks further away from half-filling and 
hopefully does nothing interesting. The FS of the growing species will now come close to 
nesting as the field is raised. This species of polarized electrons constitutes the spinless 
fermions of the calculation, while the applied field, which adds a constant (Zeeman) energy 
per particle, will play the role of the changing chemical potential. Note however that as 
the magnetic field grows and pushes the FS closer to nesting, the / in the spinless version 
decreases. At one point the system should enter the inhomogeneous CDW phase, which here 
means the spin density also will oscillate with a momentum slightly different from (ir, 7r). 
Eventually the system will jump by a first order transition to the (ir, n) state. In other 
words, when the chemical potential (i.e.e, applied magnetic field) 

falls within a window, the system density locks at half-filling, taking particles from the 
reservoir, which in this case is the other species of spins. We should see a sudden increase 
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in polarization which should then remain fixed as the field is increased, since the up and 
down densities remain locked. (Unlike in the case of an ideal reservoir which can donate 
any number of particles at one energy, there will be a cost of energy when a macroscopic 
number of fermions have to be converted from the shrinking FS to the growing one. This 
will decrease the region of stability of the commensurate (ir, n) CDW phase.) If the field is 
raised further, we are effectively changing the sign of / and so the same picture will appear 
in reverse: a first order transition to the inhomogeneous state followed by a second order 
transition to the Fermi liquid. 

By design, the grand canonical picture applies to the above experiment, where the system 
of polarized electrons indeed is in contact with a reservoir (of opposite spin electrons). The 
free parameter is the chemical potential or the applied magnetic field, and the particle 
density is chosen by the system to minimize its energy. 

Now the ideas discussed above also apply to other problems where the number density 
is the independent variable, say when a system of regular electrons at half-filling, (as in a 
Hubbard model), is doped. Whereas if we took a commensurate system and changed its 
chemical potential, it will initially ignore it and hang on its special density, here we want 
to forcibly change the number density and see what happens. We illustrate what is to be 
done in such a case by re-expressing the above spinless fermion analysis in terms of number 
density. In other words we ask- Suppose we took a half-filled system of spinless fermions 
and doped it, what will it do? Figure 4, which is a schematic, helps us analyze this. Along 
the x-axis we vary /, the chemical potential, and cover the three phases. (The CDW phase 
refers to the state). We find (in our calculation) that it has half-filling for all / up 

to the first order transition. So we plot the corresponding n as a dark horizontal line ( 
at 1/2) going up to the transition. A slight increase in I causes a first order transition to 
the inhomogeneous state. It has a lower density n q . Further increase in I causes a decline 
in n till we reach the second order transition to the Fermi liquid. There is no jump in n. 
Hereafter n and / are related as in a free Fermi liquid: n — \ — I. (The coefficient of /, 
the density of states, is set to unity). Now we can answer the question we posed. To find 
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out what happens at any given density, choose a value for n, move horizontally to the dark 
line and come right down to the phase diagram on the / axis. Thus when n — 1/2 the 
system is in the (ir, ir) state. As n is lowered, ever so slightly, the system is stuck at the first 
order transition and there will be coexistence of the (ir, ir) CDW and inhomogeneous states, 
the ratio being determined by n. When n is lowered down to n q , the system will become 
all inhomogeneous. Thereafter we move smoothly along the dark line to the Fermi liquid 
transition and Fermi liquid phase. Thus, although a slight change in chemical potential does 
nothing to the (ir, ir) state, the slightest doping causes phase coexistence. 

The above discussion has focused on breaking the nesting symmetry by changing the 
chemical potential. This is the choice with the closest analogy to the superconductor. We 
have also studied a case where nesting is destroyed by adding a second neighbor term which 
causes additional wiggles on the FS without changing its volume. Here we found that the 
window for the inhomogeneous state shrinks. In the large r limit, it goes to zero: the second 
order transition at the origin in A space occurs when the first order transition does. It is still 
interesting that the first order transition far from the origin will be accompanied by huge 
fluctuations (at an incommensurate momentum) close to the origin. The reason the uniform 
state is more stable in the presence of this kind of modification of the FS is clear: whereas 
changing / ruins nesting at (ir, ir) uniformly over the FS, the additional term due to second 
neighbor hopping vanishes at some points on the old FS. Indeed any modification to the 
dispersion relation will necessarily be periodic and vanish somewhere, unless the periodic 
function has a constant component, as in the case of the 

chemical potential. Only in the latter case, when the new free-electron FS does not 
intersect the old FS anywhere, do we have a big window for the inhomogeneous state. 

A problem worth studying is one where both a second neighbor term and chemical 
potential destroy the nested Fermi surface, with the latter as the control parameter. 
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IV. THE RENORMALIZATION PROGRAM 



As mentioned earlier, the simple minded and intuitively appealing procedure of eliminat- 
ing thin slices on either side of the FS works in the case of problems where the FS has special 
symmetries. The reason, as explained in detail in reference 0, is as follows. When we com- 
pute a one loop graph for the /3-function, the two lines in the loop are either particle-hole (in 
the CDW case) or particle-particle (in the BCS case). To get a nonzero contribution to the 
flow, it is necessary that both lie in the thin shells being eliminated and also obey momentum 
conservation. In the BCS case if the overall momentum is zero, these lines have equal and 
opposite momenta. Given time-reversal symmetry, they have equal energies, so that if one 
lies in one or the other shell (above or below the FS) that is being eliminated, so does the 
other, no matter what the direction of the line momentum. In the CDW case, the momenta 
of the particle and hole lines differ by Q = (it, it). Since this reverses the energy for the 
nested problem, again if one lies (anywhere) in the shell below or above the FS that is being 
eliminated, the other lies in the shell of opposite energy, also being eliminated. 

In both cases, no matter how many shells we eliminate, the couplings keep flowing, with 
each factor of s reduction in cut-off giving a contribution proportional to Ins. The flow is 
off to the condensed state. 

All this changes if we break the symmetry. Consider the CDW example. Now in order 
for one of the momenta to lie in a shell of thickness dA and scatter into another one also 
being eliminated, its momentum must have a very specific direction that lies within a narrow 
range of size equal to the shell thickness divided by bandwidth. Thus in the shell elimination 
game, there will be no flow. This is not an insurmountable problem. In reference [|IJ another 
scheme called the Field Theory Scheme (FTS) is invoked. Here one computes some physical 
quantity in the cut-off theory and sets the derivative with respect to the cut-off to zero to 
obtain the (5- function. We get the same flow as before for A >> /. The flow does change 
character as A approaches /, as it should, and indicates a transition. But we are not very 
happy with this method since the FTS is not generally to be used when any of the other 
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energies in the problem comes close to A. This is because in the FTS one tries to get 
away with just the quartic interaction and this requires that the cut-off dominate all other 
energies. In other words the smallness of ratios like p/A or q/A is what allows one to neglect 
higher operators. In the full fledged Wilson-Kadanoff scheme there is, of course, no such 
restriction, provided one works with operators of arbitrary complexity, a prospect we do not 
want to entertain. We are currently engaged in solving this problem. 



Another approach is to follow Hertz fl3fl . Here one couples the fermions to a bosonic 



variable A as in a Hubbard-Stratanovich transformation, integrates the fermions out and 
works with the effective action for A. As mentioned in Reference [0J, this is generally going 
to cause trouble since integrating out gapless degrees of freedom (anathema to the standard 
RG) can and typically will, lead to a singular action for the remaining fields. If one did this 
for the usual BCS problem, one finds that the action cannot be expanded in powers of A. 
Even with the / term present, the action is non polynomial in A since the chosen momentum 
at the transition to the inhomogeneous phase is a singular point of the susceptibility. (In 
three dimensions this is not so because of phase space. In one dimension the singularity 
is even more pronounced.) Even if the action is polynomial, it may not have an analytic 
expansion in u or q. Hertz has argued that in some cases we can still find the right scaling 
when this happens. 

Returning to the general problem, one option is to keep both the fermions and bosons 
together and not eliminate only non-singular modes. This is what one does in problems 
where gauge fields couple to fermions. In such cases the following question arises. For 
bosons the energy is measured from the origin in momentum space while for fermions it 
is measured from the FS. How are high and low energy modes to be defined? A boson 
that imparts a momentum parallel to the fermions momentum can take it out of the cut-off 
(and is a high energy boson) while the same boson, if it attaches itself to a fermion with a 
perpendicular momentum, will move it along the FS and be a low energy boson. A common 
solution is to pick a point on the fermion's FS and treat bosonic momenta in radial and 
angular directions distinctly. We find this approach to a rotationally invariant problem 
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unsatisfactory , even though for reasons we do not fully understand, it might work. On the 



other hand a recent paper by Altshuler, Ioffe and Millis |H] describes a class of problems 
where such a choice of coordinates is warranted. These authors consider a FS which is not 
rotationally invariant and only two small segments of which (parallel under translation by 
a vector Q) are important. In this case there is a preferred fermion direction and using the 
fermion dispersion relation near these points to decide which bosons are high energy and 
which are low energy seems legitimate. This is one very promising possibility we intend to 
fully digest when we resume our RG program. However, for isotropic problems involving 
bosons and fermions, we believe there is still need for improvement. In the meantime we 
can take the attitude that since no problem is really isotropic, one may begin by attacking 
the gauge problems by starting with, say an elliptical FS and trying the scheme of Altshuler 
et at 



V. CONCLUSION 

We began with systems that have quantum phase transitions at zero coupling, which 
exist due to special symmetries of the FS. Two examples given in Reference Jl| and discussed 
here were the planar superconductor (with a time-reversal invariant spectrum) and planar 
spinless fermions at half-filling (with nesting symmetry of the FS). The goal was to analyze 
systems with transitions at nonzero, but small, coupling so that a perturbative RG would 
be possible. To this end we had to come up with systems with such transitions, study them 
in mean-field to map out the phase structure, and then apply the RG to the transitions. 
The first step was accomplished by destroying the symmetries of the FS by adding a term 
/ to the hamiltonian. In the superconductor this was a parallel magnetic field and in the 
spinless fermion case it was a chemical potential, which we saw could be accomplished once 
again by a parallel magnetic field. The, second step, which led to the results emphasized 
here, is the 

mean-field analysis which showed us the following : 
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The two systems have a very similar three phase structure and phase transitions be- 
tween them. They have isomorphic formulas for energy and susceptibilities. They 
naturally have similar windows for the inhomogeneous state. Note that we do not 
expect the systems to behave identically when fluctuations are considered, since they 
break different kinds of symmetries. 

The phase at large / (or small coupling) is a Fermi liquid, which 

is polarized in the superconducting case. The next phase as / is reduced is has a 
condensate that has a nonzero momentum with respect to the canonical value, this 
being zero for the superconductor and (n, tt) for the CDW. For this reason we call 
it the inhomogeneous phase. It has gapless excitations along with the nonzero order 
parameter. We have computed some of its properties for the superconductor in Section 
III. The last phase (at small /) is a state with condensate at the canonical momentum, 
which is the BCS state for the superconductor and the Q = (tt, tt) state for the CDW. 

The transition from the Fermi liquid to the inhomogeneous state is second order in 
mean-field but with nonclassical exponents: (5 = 2 for the superconductor. The non- 
classical exponents arise because the chosen momentum for the condensate is a singular 
point of the Fermi liquid susceptibility. The transition from the intermediate state to 
the homogeneous state is first order. 

The large window for the inhomogeneous state is due to two-dimensionality. Indeed 
anisotropic systems have larger windows since they are more one-dimensional. Recall 
that for the CDW case 

^ < / < rA„ 

where A is the order parameter before any field is applied and r is the anisotropy. 
For the isotropic superconductor, set r = 1 above. 

We expect a similar enlargement of the window for anisotropic superconductors. 
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• In terms of experimental realizations, we require clean, truly two-dimensional, low Z 
systems. 

In the superconducting case we need a B (in Tesla ) of order T c 
(in Kelvin) assuming a free electron g factor. 

In the CDW case, in an idealized world where spinless fermions exist, once again 
Zeeman energies of the order 

of the gap will be needed. 

In practice, to employ the trick discussed here, stronger fields may will be needed to 

separate real electrons into two sufficiently polarized species. Here a large g factor will 
help. 

Our analysis has many limitations. First it is mean-field. The order of the transitions 
may change. (The analysis of Altshuler et al suggests both transitions may be first order.) 
However, given the large windows, the intermediate phases should not disappear due to 
fluctuations. 

Our analysis ignores impurities and spin-orbit scattering. To see the phase discussed 
here will require clean, low Z two-dimensional systems and strong magnetic fields. For the 
CDW case, where we try to move the FS using the field, a large g factor will help. It appears 
that we are at the threshold of having all these. 

Our analysis is at T = 0. This restriction was not due to any technical reasons and can 
be readily overcome. 

While our RG program has naturally led us to study this family of systems, the systems 
themselves are not new. As we mentioned, Bulaevskii has presaged many of our results for 
the superconductor. There is no doubt a similarly huge body of literature on the CDW 
problem, both theoretical and experimental. Rather than focus on individual systems, we 
have taken a slice through them and emphasized the common thread that runs through 
many of them. The next step is to go through individual cases that have been studied 
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experimentally and see if they could belong to this family, i..e, see if they were clean enough, 
two-dimensional enough, satisfy the various assumptions made in the derivation etc. For 
example do our consideration possibly apply to granular superconductors? [13] ]TB| Do they 
apply to CDW transitions induced by application of pressure? JT7] We expect a careful 
analysis of the huge literature (of which we have touched on a small sample just to make our 
point) will be quite involved and welcome the readers' input. In the meantime we emphasize 
that the 

experimental detection and study of the 

intermediate inhomogeneous phase in any one system will not only be fascinating in its 
own right, but also in terms of the family structure emphasized here. 
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Figure Captions 



Figure I. Linear phase diagram for (a) the superconductor and (b) CDW systems as / is 
raised. In both cases, 1$ is where the Fermi liquid susceptibility for uniform condensate 
diverges. It does not correspond to a phase transition. In the CDW case I q is when the 
Fermi liquid susceptibility diverges for momentum is the x direction. This too is not a real 
transition since the system has already condensed at I p to a state with momentum in the y 
direction. All momenta are measured from (ir, 7r)) in the CDW case. 

Figure 2. The dark lines show the anisotropic, two-branched, nested Fermi surface at half- 
filling. The vector Q) connects any point on it to any other point. When the chemical 
potential is turned on, the surface shrinks to the thin line. There is no global nesting vector. 
Figure 3. The energy of the the inhomogeneous and uniform BCS states relative to the 
polarized Fermi liquid (PFL) as a function of the order parameter A, just when the BCS 
energy (labeled q = ) relative to that of the PFL. vanishes. This is the point Ip = A /v^2 
in Figure 1. Note that that at this point the inhomogeneous state has a lower energy. 
The optimal momentum is somewhat larger than I. Given the energy scales, it clear that 
under a very slight reduction of /, the BCS state will dip below the inhomogeneous state. 
Thus we may take Ip = A /^/2 as the transition point between the uniform state and the 
inhomogeneous state. 

Figure 4. Phase structure as a function of number density n as well as chemical potential 
/. As / is raised, the system evolves along the x axis from the (ir, it) CDW state to the 
inhomogeneous state to the PFL. As a function of n, slightest doping from half-filling leads 
to phase coexistence till n comes down to n q . Thereafter the behavior is quite smooth. 
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